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Abstract. We propose a nonparametric method to determine the sign of 7 - the 
redshift evolution index of dark energy. This is important for distinguishing between 
positive energy models, a cosmological constant, and what is generally called ghost 
models. Our method is based on geometrical properties and is more tolerant to 
uncertainties of other cosmological parameters than fitting methods in what concerns 
the sign of 7. The same parameterization can also be used for determining 7 and 
its redshift dependence by fitting. We apply this method to SNLS supernovae and 
to gold sample of re-analyzed supernovae data from Riess et aZ.2004. Both datasets 
show strong indication of a negative 7. If this result is confirmed by more extended 
and precise data, many of dark energy models, including simple cosmological constant, 
standard quintessence models without interaction between quintessence scalar field(s) 
and matter, and scaling models are ruled out. 
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Recent observations of Supernovae (SNe), Cosmic Microwave Background (CMB), 
and Large Scale Structures (LSS) indicate that the dominant content of the Universe is a 
mysterious energy with an equation of state very close to Einstein cosmological constant. 
The equation of state is defined by w, the ratio of pressure p to density p, w = P/ p. For 
a cosmological constant w = —1. The observed mean value of w for dark energy is very 
close to —1. Some of the most recent estimations of w are: From combination of 3-year 
WMAP and SuperNova Legacy Survey (SNLS), w = — 0.97^q9°^[2J; from combination 
of 3-year WMAP, large scale structure and supernova data, w = — l.OGlo^oogH; from 
combination of CMAGIC supernovae analysis and baryon acoustic peak in SDSS galaxy 
clustering statistics at 2; = 0.35, w = —1.21~^q^I^; and finally from baryon acoustic peak 
alone w = — 0.8±0.18. It is evident that with inclusion of one or two sigma uncertainties 
to measured mean values, the range of possible values for w runs across the critical value 
of —1. Moreover, in all these measurements the value of w depends on other cosmological 
parameters and their uncertainties in a complex way. Reconstruction methods for 
determining cosmological parameters from observations - see [5] and references therein 
for a review of methods - usually use fitting of continuous parameters on the data 
and determine a range for w. On the other hand, the sign of w -|- 1 is more crucial for 
distinguishing between various dark energy models than its absolute value. For instance, 
if w -|- 1 < 0, quintessence models with conventional kinetic energy and potential are 
ruled out because in these models this quantity is always positive. Decay of dark matter 
to dark energy or an interaction between these components can lead to an effective 
w -|- 1 < without violating null energy condition [8]. Redshift dependence of w has 
been invoked in many of recent works. Although the quality of available data does not 
yet permit to obtain a precise measurement of variability of w, recent attempts show 
that at low redshifts this can not be very large [7], see also the results of present work 
below. 

Without loss generality the density of the Universe at redshift z can be written in 
the following form: 

^ = nmil + Zf + + ZY + fi,e(l + (1) 

Po 

where p{z) and po are total density at redshift z and in local Universe, respectively; 
flm, ^h, and fide are respectively cold and hot matter, and dark energy fraction in 
the total density at z = 0. We consider a flat universe in accordance with recent 
observations |[2J. At low redshifts, the contribution of the CMB to the total mass of 
the Universe is negligible. The contribution of neutrinos is Qi^h"^ = ^m;^/92.8 eV, 
h = ifo/lOO/cm Mpc~^sec~^, where Hq is present Hubble constant. The upper limit on 
the sum of the mass of neutrinos from 3-year WMAP is < 0.62 (95% confidence 

level) pZI]. Even if one of the neutrinos has very small mass and behaves as a warm 
dark matter, for z < 1 neutrino contribution to the total mass is < 4%, smaller than 
the uncertainty on dark matter contribution. Thus, the approximation Qm + fide ~ 1 is 
justified. 

Here we propose a nonparametric method specially suitable for estimating the sign 
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of 7. When the quahty of data is adequate, the quantity A{z) defined in can also be 
used to fit the data and to measure the value of 7. The expression nonparametric here 
is borrowed from signal processing literature where it means testing a null hypothesis 
against an alternative hypothesis by using a discrete condition such as jump, sign 
changing, etc., in contrast to constraining a continuous parameter (see e.g. ^). We 
show that geometrical properties of A{z) are related to the sign of 7 and we can detect 
sign without fitting a continuous parameter. 
Assuming a constant w 

-f = w + l (2) 

and it can be easily shown that in this case: 

A{z) ^ ^ ^ -n^ = 7l^.e(l + zf'-^-'^ (3) 
3(1 + zj'^po dz 

If the equation of state of dark energy varies with redshift, we can assume that 7 in 
([1]) depends on redshift. In this definition 7 is related to the equation of state of dark 
energy according to: 

liz) = r-- / dz (4) 

' ln(l + ^) l + z' ^ ^ 

and: 

A{z) = n,, {j{z) + il + z) ln(l + z)^) (1 + zf^''-''^ (5) 

When w is constant (jl]) and ([5]) approach respectively to ([2]) and ([3]). If w{z) is a slowly 
varying function of redshift, we can parameterize it as: 

w{z) = Wq + WiZ (6) 

where Wq and Wi are constant, and \wi/wo\ < 1. This leads to: 

-f{z) = l+Wo + Wi{z/\n{l + z) -1) (7) 

By expanding the logarithm term at low redshifts it is straightforward to see that at 
lowest order 7 ^ 1 + wo + wiz/2. Therefore 7(z) varies even more slowly than w. 
Using ([7]) and ([5]), we can find the explicit expression of A{z) as a function of constant 
parameters: 

A(z) = n.(7(.) + """'7J(lt.!"''' )<^ + 



^l],e(7W+^i(y + ...)j(l + ^r^-'^ (8) 

This shows that at low redshifts the contribution of derivative term is quadratic in z and 
negligible in comparison with 7(2;). Therefore, at low redshifts the sign of A{z) follows 
the sign of 7(2) even for a redshift-dependent equation of state. On the other hand, 
for z < 1 this parameterization is equivalent to models IV and IV of Serra et al. [7j, 
which are barely constrained and are consistent with zero or small Wi. This confirms 
that the assumption of a constant or slowly varying 7 at low redshifts in ([3]) is a valid 
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approximation for our Universe. Therefore in this work from now on we consider only 
this case, unless the redshift dependence of 7 is explicitly mentioned. 

Similar expressions can be obtained for non-standard cosmologies such as DGP[T2] 
model and other string/brane inspired cosmologies [13j. It is also possible to find 
an expression similar to ([3]) for non-flat FLRW models and without neglecting hot 
matter. The left hand side would however depend on flk and flh, and therefore would 
be more complex. Nonetheless, when their contribution at low redshifts are much smaller 
than cold matter and dark energy, the general behaviour of A{z) will be the same as 
approximate case studied here. 

In summary, the right hand side of ([3]) has the same sign as 7 if 17(2;) | 3> 
(1 + z) ln(l + 2;)|7'(z)|. Moreover, the sign of A{z) derivative is opposite to the sign 
of 7 because due to the smallness of observed 7, the term 7 — 1 is negative (if 7 depends 
on z this is true only for low redshifts). This means that A{z) is a concave or convex 
function of redshift, respectively for positive or negative 7. In the case of a cosmological 
constant A{z) = for all redshifts. This second feature of expression ([3]) is interesting 
because if Qm is not correctly estimated, A{z) will be shifted by a constant, but its 
geometrical properties will be preserved. 

The left hand side of expression ([3]) can be directly estimated from observations. 
More specifically, flm is determined from conjunction of CMB, LSS, and supernova type 
la data, and at present it is believed to be known with a precision of ~ 5%. At low 
redshifts, the derivative of the density is best estimated from SN type la observations. In 
the case of FLRW cosmologies, the density and its derivative can be related to luminosity 
distance Di, and its first and second derivatives: 

3(1 + zypo^dz i(f^-^)3 

A = (l + .)i/oy^'^ , H^{z) = ^-^p{z) (10) 

It is remarkable that the right hand side of (Q depends only on one cosmological 
parameter, Hq. Moreover, similar to an uncertainty on Qrn, Hq uncertainty rescales 
B{z) similarly at all redshifts, but does not modify geometrical properties B{z) and 
A{z). Di can be directly obtained from observed luminosity of standard candles such 
as supernovae type la. In the case of LSS observations where the measured quantity 
is the evolution of p with redshift, B{z) is measured directly up to an overall scaling 
by pq. This does not change the geometrical properties of B{z) and A{z) either. The 
scaling by a positive constant does not fiip a convex curve to concave or vis-versa. 
In summary, uncertainties of Qm and Hq do not affect the detection of the sign of 7 
through geometrical properties of A{z). This is quite different from fitting methods. 
They are sensitive to all numerical parameters Hq, Q^, ^de, and w in a complex way, 
usually through a non-linear equation such as chi-square or likelihood. This makes the 
assessment of the infiuence of uncertainties on parameter estimation, and specifically on 
the determination of the sign of 7 very difficult. 
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When this method is apphed to standard candle data such as supernovae type 
la, one has to calculate first and second derivatives of Di. Numerical calculation of 
derivatives is not however trivial. To have a stable and enough precise result, not only 
the data must have high resolution and small scatter, but also it is necessary to smooth 
it. Therfore, to test the stability of numerical calculations and the method in general, 
we also apply it to simulated data. The details of numerical methods is discussed in the 
Appendix. 

Examples of reconstruction of simulated data are shown in FigUl It is evident 
that despite deformation of the reconstructed curves due to numerical errors and 
uncertainties, specially close to lower and upper redshift limits, the difference between 
convexity of curves for positive and negative 7 is mostly preserved and can be used, 
either by visual inspection or by using a slope detection algorithm, to find the sign of 
7. The simulated data is however much more uniform - has less scatter - than presently 
available data, see FiglH Therefore the reconstruction from real data is more prone 
to artifacts than simulated ones. FigH] shows that although at mid-range redshifts the 
data follow a curve similar to what is expected from FLRW cosmologies, deviations 
appear close to boundaries and at high redshifts where the quahty of data is worse. We 
attribute these to the artifacts of reconstructions. 

Visual inspections or slope detection lack a quantitative estimation of uncertainties 
of the measured sign for 7. In signal processing usually binomial estimation of the 
probability or optimization of detection[10] is used to assess uncertainties. However, 
in most practically interesting contexts in signal processing, the signal is constant and 
uncertainties are due to noisy. In the cosmological sign detection discussed here the 
observable A{z) is both noisy and varies with redshift. Therefore, binomial probability 
and similar methods are not applicable. For this reason we try another strategy, specially 
suitable for this cosmological sign detection task. 

The null hypothesis for dark energy is 7 = 0. Assuming a Gaussian distribution for 
the uncertainty of reconstructed A{z) from data and from simulated data with 7 = 0, 
for each data-point we calculate the probability that the data-point belongs to the null 
hypothesis. To include the uncertainty of data, we integrate the uncertainty distribution 
1-sigma around the mean value: 



where Ai and cTj belong to the i*^ data-point, and A^^ and (Xoi belong to simulated null 
hypothesis model at the same redshift. Averaging over Pi gives P, an overall probability 
that the dataset corresponds to the null hypothesis. As 7 = is the limit case for 7 > 0, 
P is also the maximum probability of 7 > 0. 

We have applied this algorithm to two supernova datasets: published data from 
Super Nova Legacy Survey (SNLS)[5], and to supernovae with z < 0.45 from gold sample 
of re-analyzed data by Riess et a/.2004[lj. The reason for using only the low redshift 
subset of the latter compilation is that scatter and uncertainties of the peak magnitude 
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at higher redshifts is too large and reconstructed A{z) is instable. FigH] shows A{z) 
obtained from these data as well as from simulated SNe as described in the caption. 
Simulated standard sources are put at the same redshifts as real data to make simulated 
samples as similar to data as possible. In both datasets the probability of 7 < or 
equivalently w < —1 is larger than 70%. The SNLS data is consistent with a 7 as 
small as ~ — O.4J. There is however significant deviation from a smooth distribution for 
z < 0.1 and z > 0.5. We attribute this to large scatter of the data at these redshifts see 
Figl2]-a. The fact that at intermediate redshifts in these datasets A{z) follows closely 
models with constant 7 is a demonstration of very small variation of 7 in these redshifts, 
consistent with our arguments above. 
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Figure 1. A{z) from 117 SNe of the SNLS (a-purple curve), and from 88 SNe with 
z < 0.45 recompiled and re-analyzed by Riess et aZ.2004[T] (b-purple curve). Error 
bars are 1-sigma uncertainty. In (a), green, orange, yellow, and light green curves are 
the reconstruction of A{z) from simulations for 7 = —0.2, —0.06,0.6,0.2, respectively. 
The probability of null hypothesis (7 = 0) is P = 0.27, therefore the probability of 
7 < 0, 1 — P = 0.73. Light grey and cyan curves are theoretical calculation including 
the uncertainty on fide, respectively for 7 = ±0.06, ±0.2. In (b), the pink curve is 
simulated 7 = model. For this dataset 1 — P = 0.75. The dark grey straight line is 
the theoretical expectation for A{z) when dark energy is a cosmological constant. For 
all models Hq = 73 km Mpc~^ sec~^ and Qde = 0.77. 



To see whether the large negative 7 concluded from the SNLS data is due to 
scattering and/or reconstruction artifacts, we also apply the same formalism to a subset 
of these data with z < 0.45. The result is shown in Figj2]-b along with simulations in 
the same redshift range. The bump at very low redshifts in Figil] does not exist in this 
plot, and therefore we conclude that it has been an artifact of numerical reconstruction. 
Although A{z) distribution in this dataset is also convex and the probability of 7 < 
is larger than 90%, it does not have the same slope as any of constant w models. It 
seems that low and high redshift sections of the curve correspond to different values of 
7. For z < 0.15, A{z) is close to theoretical and simulated data with 7 = —0.2. For 
z > 0.25, A{z) approaches the values for larger and even positive 7. Such behaviour 

J Here we concentrate on the sign of 7 and don't perform any fitting to obtain its value. Estimations 
are from comparison with simulations. 
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does not appear in FigUl Giving the fact that the number of available data points with 
z > 0.25 in this subset is small, the most plausible explanation is simply numerical 
artifacts. Alternative explanations are the evolution of 7 with redshift or the use of 
under-estimated value for Qrn- If the latter case is ture, the value of 7 must be even 
smaller than —0.2. Interestingly, the deviation from a constant w model in the former 
case is consistent, up to uncertainties, with the best estimations of the evolution of w in 
models IV and VI of Ref . [7] . This confirms the consistency of two reconstruction and 
parameter estimation methods. However, as explain above, present method shows that 
this deviation can be also due to uncertainties of Qm and not the evolution of w. The 
simplicity of dependence on the cosmological parameters in this method permits to see 
their effects more explicitly than in fits. With a data gap in 0.15 < z < 0.25 interval 
and only 58 supernovae in this subset, it is not possible to make any definite conclusion 
about the behaviour of this data. We should also mention that SNLS supernovae with 
z < 0.25, and at 0.25 < z < 0.4, and z > 0.4 have not been treated in the same way[5]. 
It is therefore possible that some of the observed features are purely artifacts of the 
analysis of the raw observations. 

FigHJ-b shows A{z) from gold SNe sample recompiled by Riess, et al.\l]. It is also 
consistent with 7 < with a probability ~ 75% at 1-sigma and ~ 66% at 2-sigma. 
This plot shows that the value of 7 estimated from this data is ~ —0.06, larger than 
estimation from SNLS data. 
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Figure 2. a: /if, from SNLS supernovae, data (magenta), smoothed distribution (blue). 
Although this distribution look quite smooth, even small sudden variations can lead to 
large variations in derivatives, b: A{z) for SNLS supernovae with z < 0.45. Definition 
of curves is the same as in FigH) For this dataset 1 — P = 0.93 when cosmological 
parameters are the same as Fig[T] If fide — 0.73 is used, 1 — P = 0.96. 



The reason for the difference between estimated values for 7 from SNLS and Riess, 
et a/, compilation here is not clear because value of w reported in Ref. [1] and Ref [5] 
are consistent. Nonetheless, the estimated values are in the range reported by other 
works [2] [3]. Most probably this difference is related to different scatter and uncertainty 
of these datasets, and the fact that low and high redshift data are not treated in the 
same way. Recent claims about contamination of supernova type Ic and the effect of 
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asymmetric explosion in the lightcurve of supernova type Ia[Tlj, and possible differences 
between low and high redshift super novae [15] can not explain the differences either, 
because they affect both surveys in the same way. Despite these discrepancies, both 
datasets are in good agreement about negative sign of 7. 

These results highlight various shortcomings in both datasets used here. Our first 
remark is the large gap in redshift distribution of observed supernovae in redshift range 
0-1 ~ -2 ^ 0-3. Both datasets have less than 6 supernovae in this range and nothing 
in 0.15 < 2 < 0.25. This is not an important issue for finding redshift distribution 
of Di, but redshift gap becomes very important when derivatives of Di are calculated. 
The lack or rareness of supernovae data in this redshift range is related partly to the 
history of star formation in galaxies[T6], and partly to optimization of surveys p!7] for 
detection of very low or high redshift supernovae which decreases the probability of 
detection of mid-range SNe. Sloan Supernova Survey [1 8 j is optimized to detect SNe in 
0-1 ^ -2^ ~ 0-35 and should provide the missing data in near future. 

Our second observation is a large scatter in both datasets around redshift ~ 0.55, 
see Figj2]-a. This leads to a large scatter in the numerical determination of derivatives 
in and makes the results unusable. In future observations the reason of this large 
scatter should be understood. 

From theoretical calculation and simulations shown in FigiU one can also conclude 
that with present uncertainties of cosmological parameters, the most important redshift 
range for determining the equation of state of dark energy is z < 0.8. Higher 
redshift supernovae are only interesting if w varies significantly with redshift. Although 
technical challenges, understanding of physics of supernovae and their evolution [T5j. 
and applications for other astronomical ends make the search of supernovae at larger 
distances interesting, they will not be in much use for determining the equation of state 
of dark energy, at least not at the lowest level which is the determination of redshift 
independent component of w. 

On the other hand, improvements in numerical techniques and algorithms would 
lead to better measurements. One of the possibilities in this direction is the application 
of an adaptive smoothing algorithm with variable degrees of smoothing depending on 
the amount of scatter in the data. More sophisticated smoothing algorithms also have 
been proposedpjj. We postpone the apphcation of these advanced methods to future, 
when larger datasets become available. 

In summary, we have proposed a nonparametric formalism to investigate the sign 
of 7 in the equation of state of dark energy. When data with better quality become 
available, fitting can be added to this method to find also the value of 7. The advantage 
of this method is its geometrical nature when the sign of 7 is searched. Moreover, its 
simple dependence on cosmological parameters makes the assessment of uncertainties of 
the value of 7 easier. This method is suitable for applying to observations in which the 
total density of the Universe and its variation are directly measured such as standard 
candles and galaxy surveys, but not to integrated observables such as CMB anisotropy. 
By applying this method to two of largest publicly available supernovae datasets we 
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showed that they are consistent with a w < —1. Present data is not however enough 
precise to permit the estimation of Iti? + 1| with good certainty. With on-going projects 
such as SNLS, Supernova Cosmology Project ^2Qj, and SDSS SNe survey, and future 
projects such as SNAP and DUNE, enough precise datasets should be available soon. 

Appendix: For standard candles, A{z) must be calculated from luminosity 
distance which is related to mag nitude: A/^o = 10''''/^ Dq is the dist ance for which 
the common luminosity of standards is determined and depends on Hq. Therefore, fif, 
needs a correction if a different Hq is used[5]. For simulated data, Di is calculated from 
(fTOj) and fif, is determined from definition above, then a random noise with a standard 
deviation of 3% is added to the magnitude. 

Expression (Q for B{z) contains the first and the second derivatives of Di which 
must be calculated numerically from data. It is however well known that a direct 
determination of derivatives leads to large and unacceptable deviation from exact values. 
One of the most popular alternatives is fitting a polynomial around each data point and 
then calculating an analytical derivative using the polynomial approximation in place 
of the data. We use this approach to determine derivatives of /if, and Di. In addition, 
before applying this approach, we smooth the distribution of magnitudes using again 
the same polynomial fitting algorithm. In FLRW cosmologies the redshift evolution 
of the luminosity distance is very smooth. Therefore, a second order polynomial for 
smoothing is adequate. Fitting is based on a right-left symmetric, least algorithm, 
and for this we have implemented Ifit function of Numerical Recipes^21j. By trial and 
error we find that 19 data-point fitting gives the most optimal results regarding the 
number of available data points and their scatter. Close to boundaries however less 
data point for fitting is available in one side of each point, and therefore the fitting 
is less precise. The artifacts discussed above are mostly related to this imprecision of 
calculation near boundaries. In the present work no adaptive smoothing is applied. In 
addition to smoothed data and their derivatives, the function Ifit calculates a covariant 
matrix for uncertainties of parameters (derivatives). We use diagonal elements as 1- 
sigma uncertainty of the smoothed data and its derivatives. The uncertainty of A{z) is 
calculated from the uncertainty of terms in ([3]) and (Q using error propagation relation: 
For /(xi,X2, . . .), o'j = J^i'^xii^f /^^i)^- Smoothed terms and parameters in A{z) are 
considered as independent variables with their own uncertainty. 
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